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Galileon inflation is a radiatively stable higher derivative model of inflation. The model is deter- 
mined by a finite number of relevant operators which are protected by a covariant generalization of 
the Galileon shift symmetry. We show that the nongaussianity of the primordial density perturba- 
tion generated during an epoch of Galileon inflation is a particularly powerful observational probe 

of these models and that, when the speed of sound is small, /nl can be larger than the usual result 
* — a 
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I. INTRODUCTION 



Inflation is an era during which the cosmological scale 
factor a(t) satisfies d 2 a/dt 2 > as a function of cos- 
mic time t, growing by a factor e^. It is a familiar 
idea that a successful implementation of inflation, by 
which we mean obtaining sufficiently large N, requires 
the inflationary Lagrangian density C to have an approx- 
imate shift symmetry — an invariance under the transla- 
tion <j) — > S c (j) = <f> + c, where <p is the inflaton field and c 
is a constant. 

Underlying this idea is the observation that slow-roll 
implies the shift symmetry is broken only mildly, both in 
the action and the equations of motion. In terms of the 
inflationary potential, V(</>), and the Planck mass, Mp, 
these slow-roll conditions are typically expressed using 
the parameters e and 77, which satisfy 2e = Mp(V' /V) 2 
and 77 ee Mp/" jV. 

How large a breaking can be acceptable? We expect 
any effective description to be valid only for a field excur- 
sion at most of order Mp, before renormalization group 
flow introduces new physics which changes the descrip- 
tion. Requiring the first- and second-order fractional 
variation in the potential, V (</>), to be small over an ex- 
cursion of this magnitude yields 

5 c ]xlV ~ v / e< 1, and S 2 In V ~ 77 - e < 1, (1) 

where we have indicated approximate relations in terms 
of the slow-roll quantities e and 77. Therefore slow-roll 



inflation, defined by the conditions e <C 1 and \r/\ <C 
1, entails tuning V(<p) so that only very mild breaking 
occurs even over large variations in field value. 

This tuning has an important consequence. Creminelli 
observed that, once V{<f) has been adjusted to satisfy 
Eq. |T]) , we can add any operator invariant under the shift 
symmetry without spoiling the property of successful in- 
flation There is a large class of such operators, con- 
structed by applying any combination of derivatives to 
the inflaton field <f>. This yields V</>, VV</> and higher gra- 
dients, with indices contracted in arbitrary combinations. 
It follows that the most general local, diffeomorphism- 
invariant action for coupled to Einstein gravity which 
is invariant under the shift symmetry can be written Q 
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— ^i? + £ M (V0, VV0, . . .) 
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In this paper, we study how inflation can be realized 
in theories of the form ([2])- In comparison with slow- 
roll inflation using canonical kinetic terms there are new 
difficulties, associated with the appearance of unstable 
'ghost' states and stability under radiative corrections. 
One class of ghost-free, radiatively stable models with an 
interesting inflationary phenomenology has been widely 
studied. These are the Dirac-Born-Infeld (DBI) mod- 
els [HI]- However, another class of ghost-free models 
has recently been constructed 0, 0] , based on a so-called 
'Galilean' symmetry to be defined in <|TT] below. This 
Galileon field has been shown to arise naturally in theo- 
ries of massive gravity without ghosts in their decoupling 
limit @| . 

Galileon models describe an effective s/iort-distance 
theory associated with a modification of gravity on large 
scales, and have mostly received attention as models of 
dark energy. However, it is equally possible that they 
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can be used to describe cosmological evolution, and have 
been investigated for this purpose by numerous authors 
0. Ali et al. obtained constraints on the parameters of 
the model (to be discussed in ijlVp using the evolution 
of perturbations during matter domination [10]. Later, 
Creminelli et al. used an effective theory of this type 
to propose an alternative to the usual inflationary mech- 
anism for generating a primordial density perturbation 

n 

Kobayashi et al. proposed a kinetically-driven model 
of inflation supported by a Galileon-like field, which 
they called "G-inflation." In this model, Cm was taken 
to be of the form C M = P(X, (j)) + F(X, <j>)U<j>, where 
X = — V AI 0V M (/> [lH- Kobayashi et al. gave background 
solutions and studied the two-point correlation function 
of perturbations; three-point correlations in this model 
were later considered by Mizuno & Koyama |l3|. As we 
will explain, although this scenario is qualitatively similar 
to a Galileon model, the "G-inflation" picture described 
in Ref. [l2| relies on a hard breaking of the Galilean sym- 
metry. This symmetry is essential to stabilize the theory 
against quantum corrections; if it is broken, there is no 
protection against the radiative generation of higher op- 
erators. These not only spoil the predictions of the the- 
ory, but lead to ghost-like instabilities. 

In this paper we take a different approach. We argue 
that it is important to retain the Galilean symmetry, as 
far as possible. It will turn out that a non-trivial back- 
ground geometry introduces soft breaking terms, which 
lead to an interesting phenomenology. However, Galileon 
models which preserve sufficient symmetry represent an 
alternative to the DBI Lagrangian as a ghost-free, ra- 
diatively stable higher-derivative inflationary model. In 
comparison with the studies of Kobayashi et al. and 
Mizuno & Koyama we work with the most general co- 
variant completion of the Galileon-invariant Lagrangian. 

We refer to an epoch of inflation driven by a Galilean- 
invariant field (perhaps broken by a nontrivial back- 
ground) as 'Galileon inflation.' We calculate the bis- 
pectrum nongaussianity parameter /nl [3 and show 
that Galileon models can have an observational signa- 
ture which is distinct from both DBI and canonical in- 
flation. Indeed, Galileon models occupy a previously 
unexplored regime in the effective field theory of infla- 
tionary perturbations, where observable nonlinearities 
may be much larger than those obtained from even the 
highly nonlinear DBI Lagrangian. This difference arises 
because the cubic interactions are described by a com- 
bination of dimension-six and dimension-seven opera- 
tors, rather than the dimension-six operator alone, which 
would generically dominate. 

In what follows we work with signature (—,+,+,+) 
and allow Greek indices fx, v, . . . , to range over space- 
time indices, whereas Latin indices i, j, . . . , label purely 
spatial indices. We choose units in which the Dirac con- 
stant, h, and the speed of light, c, are set to unity, and 
express the gravitational coupling in terms of the reduced 
Planck mass, Mj = (8nG)~ 1 / 2 . Spacetime indices in 



complementary upper and lower positions are contracted 
with the spacetime metric, and □ = V^V^. We reserve 
d to denote a derivative with respect to purely spatial 
indices. Spatial indices are contracted using a Kronecker 
delta, so that d 2 — didi. 



II. INFLATION 

In this section we briefly review the conditions Cm 
must satisfy to constitute a well-defined, predictive de- 
scription of an inflationary phase. 

Unitarity. An arbitrary action of the form @ contains 
derivatives of second-order and higher applied to the el- 
ementary fields, and therefore leads to equations of mo- 
tion of third-order or higher. Unless an infinite number of 
derivatives are present such theories do not usually pos- 
sess a well-defined Cauchy problem [l5j ] , which manifests 
itself in the appearance of ghost states. At the quantum 
level the result is a pernicious loss of unitarity. Accord- 
ingly we should restrict our attention to special choices of 
Cm which lead to second-order equations of motion . 

Quantum fluctuations. It has recently been argued 
that the choice Cm = F(X, 4>)\^4> (used in Ref. [Hj]) leads 
to second-order equations of motion for any choice of 
F [l7j ] . and therefore unitary evolution as a quantum 
field theory. For this reason, one could contemplate such 
models as candidates for an inflationary action of the 
form ©. 

However, unitarity is not the only condition for pre- 
dictivity. A typical F can be represented as a power 
series in X and (j), with some specified coefficients. But 
we should recognize that in the absence of other sym- 
metries there is no way to protect these coefficients from 
quantum corrections. An important example is the DBI 
Lagrangian Cm 12 /((^)vl-X for an arbitrary f(<p). In 
the limit X ~ 1, fluctuations around the background 
solution have small sound speeds c s « 1 and acquire 
large nongaussianities of order /nl ~ cj 2 0, [T^]. But 
at X ~ 1 the square root formally receives contribu- 
tions from all powers of X with coefficients in precisely 
defined ratios, which in principle are susceptible to dis- 
ruption by renormalization. The prediction /nl ~ cj 2 
is believed to be trustworthy only because a symm etry 
forbids large renormalizations of these coefficients |3.ll9j. 
forcing quantum corrections to involve the £wo-derivative 
combination VV</>. A similar symmetry protects fluctu- 
ations around models that exhibit Galilean symmetries 
(to be defined below), such as those that arise in the 
decoupling limit of either massive gravity or the Dvali 
Gabadadze-Porrati (DGP) model, and can even be ex- 
tended to theories containing multiple fields [g, [D, |20h23| . 
We discuss these non-renormalization properties in mill 

The lesson we wish to draw from these examples is that 
our calculations are likely to be trustworthy only if the 
functional form of Cm is protected from large renormal- 
izations by the presence of a symmetry. In seeking to 
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generalize the DBI and Galileon-type Lagrangians, there 
are essentially only two choices. 

(a) The shift symmetry — > <f>+c is the only symmetry 
of the action. In this case we must take the leading 
relevant operators to be the quadratic mass and ki- 
netic terms, 4> 2 and (V(/>) 2 , and terms with more 
derivatives should be treated as irrelevant in the 
technical sense. Corrections to the Gaussian theory 
arise from at most a few non-renormalizable oper- 
ators, leading to very small observable departures 
from purely Gaussian fluctuations [I], [24| • 

(b) The shift symmetry is combined with a second sym- 
metry, which protects the form of any derivative 
interactions. Such a symmetry cannot commute 
with the generators of the Poincare group, because 
it necessarily mixes derivatives of different orders. 
Therefore it must be nonlinearly realized as a non- 
factorizable extension of the Poincare group. 

Choice |(a)| leads to canonical slow-roll inflation with 
very small corrections, and is well-understood. Interest- 
ing generalizations must therefore make use of |(b)| and 
can be classified according to their choice of nonfactor- 
izable extension. Only a handful of such extensions are 
known, of which one is the conformal group. Others can 
be constructed beginning with the Poincare, de Sitter or 
anti-de Sitter groups in dimension greater than four and 
applying Wigner-Inonii contractions. Following this rea- 
soning one arrives at a short list of possible Cm which are 
protected by symmetries. The choices are DBI, warped- 
DBI, Galileon and conformal- Galileon models, together 
with their higher-dimensional generalizations 0| [HJ . 
In this paper we focus on the general class of (conformal-) 
Galileon models. In such models the protecting 'Galilean' 
symmetry is an extension of the shift symmetry to in- 
clude spacetime translations, <j> — > 8 s <f> — + b^x^ + c, 
where 6 M and c are constant. 

Observables. Eq. ([2} implies a large degeneracy among 
Cm which realize given values of e and rj, and therefore 
the amplitude p 1 / 2 and spectral index n s of the primor- 
dial density perturbation (2(| . Over the last decade, sig- 
nificant effort has been invested in developing observ- 
ables which distinguish one choice of Cm from another 
[H. fl8ll24l . [27l |. As we now explain, the most powerful fam- 
ily of such observables is the amplitude of three-, four- 
and higher n-point correlations. 

How are inflationary models to be distinguished? The 
renormalization programme has taught us that predic- 
tions extracted from quantum field theories, such as 
those governing inflation, express measurable quantities 
in terms of a finite number of experimental inputs — 
sometimes expressed as "observables in terms of observ- 
ables." In canonical inflation there are two relevant op- 
erators, (V0) 2 and <fi 2 , which each require a parameter 
to be extracted from experiment. Typically these are 
(H/Mp) 2 /e and 77, extracted from V 1 / 2 and n s . The ratio 
(H/Mp) 2 alone can be determined from the tensor-scalar 



ratio r and fixes a particular de Sitter background geom- 
etry. In a noncanonical model of the form (0) , the action 
for fluctuations may exhibit spontaneous breakdown of 
Lorentz invariance, giving <fi 2 and (<9</>) 2 independent co- 
efficients. This breaking is measured by the sound speed 
c 2 . The gravitational coupling, Mp, is measured by ter- 
restrial gravitational experiments. 

We conclude that no measurement involving "P 1 / 2 , n s , 
and r alone can tell us about the operators present in 
Cm, although c s can diagnose whether <\? and (dcf>) 2 have 
independent coefficients. To distinguish one choice of 
Cm from another, we must ask whether further observ- 
ables such as the bispectrum nongaussianity, /nl, and 
the trispectrum observables tnl and c/nl [28| can be ex- 
pressed in terms of the input parameter set. Maldacena 
showed that, in a single-field model, this is true for /nl 
[29l ] unless an irrelevant operator associated with some 
higher slow-roll parameter is unexpectedly large. The 
same is true for the single-field trispectrum observables 
.'!() . It is for this reason that nongaussianities are such a 
stringent test of the single-field framework, because mea- 
surements of /nl, tnl and <?nl cannot simply be ab- 
sorbed into a new parameter of the model. 

Similar reasoning applies to noncanonical models de- 
scribed by @ . There is a simplification in theories which 
are not approximately canonical inflation, to be discussed 
in mV CI below, because the influence of gravity becomes 
negligible 0, H3, Hi- 33] . Therefore observables involve 
the parameters of Cm alone, and do not mix with Mp. 
Different choices of Cm may require a larger or smaller set 
of input parameters to express their predictions. How- 
ever, once these are measured, the remaining observables 
are uniquely determined. Therefore these express gen- 
uine differences between one Cm and another. For ex- 
ample, in DBI inflation there is a firm prediction for 
/nl once c s has been determine d, p erhaps using the 
tilt of the tensor power spectrum [3J| . We will show in 
^IV El below that there is a different prediction for /nl in 
Galileon models, making the nongaussianity of the pri- 
mordial density perturbation a particularly powerful ob- 
servational probe of Galileon inflation. The situation is 
less favourable if the form of Cm is not protected by a 
symmetry as with Cm ^ F(X, </>)□</> for general F. In 
that case, an infinite number of experimental inputs are 
required to determine F and the theory loses predictivity. 

Experimental limits. Measurements of the basic input 
parameters V 1 / 2 and n s are now quite precise [35j . and 
will improve further once data arrive from the Planck 
satellite. Limits on the bispectrum vary depending on 
the momentum dependence, or 'shape' 13611 . Present-day 
limits are moderately constraining [321 l35l. |3?| . and data 
from the Planck satellite is expected to furnish very strin- 
gent constraints. For this reason it is realistic to antici- 
pate the role of /nl as a constraint on inflationary mod- 
els. Limits on the trispectrum parameters are presently 
weak [3^ |. Smidt et al. estimate that an optimistic fu- 
ture experiment may constrain /nl to ±3, tnl to ±225 
and #nl to ±6 x 10 4 . It is not yet clear whether the 
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amplitude of higher correlations can ever be significantly 
constrained, so these parameters are likely to exhaust the 
available input parameters for at least the near future. 

If Cm requires more input parameters than can 
be measured, the model again becomes unpredictive. 
This contrasts with accelerator experiments where — in 
principle — there is typically no limit to the number of 
parameters which can be taken from experiment. In this 
paper we wish to emphasize the distinction of principle 
between a predictive theory, where a symmetry picks out 
a finite number of invariant operators whose coefficients 
fix the other observables, and unpredictive theories where 
quantum effects imply that an infinite number of coeffi- 
cients must be retained. In practice, however, a failure 
of prcdictivity for whatever cause is problematic. 



III. BREAKING THE SHIFT SYMMETRY 

An exact shift symmetry in Eq. (0) makes the value 
of 4> unphysical. To proceed we softly break the symme- 
try by hand, giving <f> a dynamical vacuum expectation 
value. We first argue that this breaking is necessary be- 
cause we would otherwise obtain 4> — 0, preventing the 
existence of a conserved curvature perturbation £ at long 
wavelengths. 

The mildest possible breaking is to introduce a poten- 
tial V = Vq — X 3 4>, where Vq is a constant and A has en- 
gineering dimension [mass]. Taking the decoupling limit 
Mp — > oo while keeping the de Sitter background fixed, 
so that 3ff 2 Mp = Vo, we recover a scalar field theory 
on exact de Sitter space. In this limit, a term linear in 
<f> does not break the shift symmetry because the action 
transforms as a total derivative, 



a(tf 



1 



-di [a(t) 3 cx l 



(3) 



We conclude that with dynamical gravity the symmetry 
must be broken softly in the sense that its effects are 
suppressed by powers of 1/Mp. 

Because the shift symmetry is exact in the decoupling 
limit, we conclude that the equation of motion for <f> can 
be written as a current conservation equation [Hj 



equations of motion.) In the unbroken case where A = 
we would obtain <f> = 0. 

At first, it might seem surprising that inflation in this 
scenario is possible only in the presence of a potential 
A ^ 0. It is now well-understood that models of "k- 
inflation" type [27| can achieve an accelerated regime in 
purely kinetic scenarios, where Cm = Cm(X). The reso- 
lution is that the equations of motion admit solutions of 
two types: a transient decaying contribution — which, in 
perturbation theory, is the decaying mode; and a domi- 
nant growing-mode solution, which is the one we consider 
above. Purely kinetic fc-inflation models make use of the 
normally transient decaying mode to source the back- 
ground expansion. In what follows we shall work with 
the more conventional solution. 

Nonrenormalization of mass. The necessity of in- 
cluding a nontrivial potential is problematic. Since po- 
tential terms are not invariant under the shift symmetry, 
they must be treated as irrelevant operators which induce 
small corrections. But operators of mass dimension less 
than four would typically receive large quantum correc- 
tions and, if present, such large renormalizations destroy 
the radiative stability of Cm , making the model of limited 
interest. In practice we cannot avoid this issue because a 
potential term of at least quadratic order is required to 
describe a graceful exit from the inflationary phase. 

Remarkably, the operators <fi and <fi 2 are protected 
by a nonrenormalization theorem and therefore can be 
treated consistently as irrelevant deformations of the 
shift-invariant Lagrangian (2lj |. This conclusion may be 
reached most straightforwardly by considering the one- 
particle irreducible effective action, T[<p], 



(6) 



expiry] = J Dtp exp jiS[0 + ip] 

We separate T[cp] into its classical part, which coincides 
with the classical action S[4>], and a quantum part T q , 
writing T[4>] =S[<p] + T q [4>]. Therefore 

expir q [0] = J Dtp exp j iS[<f> + ip] - iS[4>] 



(7) 



= j t + 3HJ t = A 3 , 



where is the Nocthcr current, 



J" 



dC 



M 



dC 



M 



(4) First, consider an action of the form 



(5) 



Eq. ([J} admits the exact solution J* = \ 3 /3H, which 
could equally have been obtained by passing to the over- 
damped limit J 1 1 J 1 — > 0. Since J* is constant, the shift 
symmetry allows us to infer that <p is also constant and 
proportional to A 3 . (To reach this conclusion we require 
the theory to have a well-posed initial value formulation, 
which is guaranteed by our restriction to second-order 



S[<t>] = S [4>] + J d 4 x ^ + -mVJ . (8) 

We assume that So[4>] is invariant under the Galilean 
symmetry, so that So^g^] = So [</>], where S g (f) — 4> + 
bfiX* 1 + c is the Galileon shift operator defined in fJTTJ It 
follows that 

expir q [0] = Jr>ip expiJ S o [0 + i/>]- S [<f>] 



f / d A x —m 2 ip 2 — ip 



6s [(f>} .sr. 



(9) 



0- 
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Therefore we conclude r q [<5 g </>] = r q [</>]. Hence, even 
though a mass term explicitly breaks the Galileon sym- 
metry, this does not induce any further operators which 
violate the symmetry at the quantum level. 

The same conclusion can be reached by analysing the 
properties of Feynman diagrams. Consider any nontrivial 
diagram with two external lines, which in principle could 
contribute to a renormalization of 4> 2 . Such a diagram in- 
cludes two field operators where the external lines attach 
to the interior of the diagram, and the structure of Cm 
implies that each of these operators carries at least one 
derivative. We conclude that the net effect of these field 
operators must introduce an overall factor of at least two 
powers of the external momenta. After expanding into 
a series of operator products, each term must contain at 
least two derivatives and cannot include the operator <j> 2 , 
contrary to our original supposition. As a trivial spe- 
cial case we recover the obvious fact that a cj) 2 operator, 
considered as a deformation of the Gaussian kinetic term 
Cm = X/2, is not renormalized. 

The same argument does not apply if wc extend @ 
to include higher powers of 4>. In particular, inclusion of 
a (j> 3 operator or higher would typically renormalize the 
coefficient of (j> 2 to ~ A 2 , where A is the cutoff of the 
theory. In such circumstances we would be obliged to 
take (j) 2 as a relevant operator, leading to strong radiative 
breaking of the supposed shift symmetry exhibited by 
Cm- Therefore inclusion of such operators is inconsistent. 

The conclusion of these arguments is that it is possible 

I 



S = 



where and G M „ and R are respectively the Einstein tensor 
and scalar curvature of the background. Typically, both 
covariantization and inclusion of a non-Minkowski back- 
ground metric will softly break the Galilean symmetry. 
It will emerge that this has important consequences for 
the phenomenology of the model. The coefficients Ci are 
dimensionless, and — as above — A is a mass scale which 
determines the naive cutoff of the theory. In practice, 
fluctuations around a nontrivial background can be valid 
up to energies somewhat larger than A if a Vainshtein ef- 
fect is operative (see Ref. |41|), discussed in more detail 
in Ref. Q. Had we included nonminimal coupling to the 
geometry in Eq. @, the curvature terms involving G M „ 
and R would have been accompanied by other geometric 
invariants, but all such terms are suppressed by powers 
of H/A (cf. Eqs. (p?U |) - (j!?5)) . where such suppressed con- 
tributions can be clearly identified). In the inflationary 



to construct an inflationary model based on a Galileon 
field, in which the shift and Galilean symmetry protects 
the form of the Lagrangian. Inflation can end, despite the 
presence of these symmetries, at least in the decoupling 
limit Mp — ► oo. Although the Galileon symmetry is bro- 
ken when coupled to gravity, the breaking terms will be 
parameterically suppressed by powers of A/Mp, which we 
assume to be small. In what follows we will allow a little 
extra freedom and work with an arbitrary potential V(4>). 
Although this will generally break the Galilean symme- 
try explicitly, the foregoing argument demonstrates that 
any operators generated in this way will be suppressed by 
powers of three or more derivatives of V(</>). We will sup- 
pose that the models used to obtain Galileon inflation all 
break the Galilean symmetry mildly in this above sense. 



IV. GALILEON INFLATION 

On a curved background, such as de Sitter space, Def- 
fayet et al. [4(J remarked that the Galileon action con- 
structed by Nicolis et al. [f| leads to unwanted higher- 
derivative equations of motion, spoiling the expected con- 
struction of a ghost-free, unitary theory. This can be 
cured using a nonminimal coupling to gravity, which Def- 
fayet et al. described as 'covariantization.' The covariant 
Galileon action, which can also be obtained from the five- 
dimensional covering theory [7j , is 



I 

regime of interest, where nonlinearities are dominated by 
the Galileon self-interactions, it will transpire that such 
terms are negligible, but for completeness we continue 
to retain the nonminmal curvature couplings required by 
covariantization. 



The action for fluctuations in the decoupling limit of 
the DGP model has 04 = 05 = 0. Constraints on the 
Ci obtained from short-distance gravitational effects were 
studied in Refs. [fl|42j for the case that Eq. flTUl) describes 
the short-distance effects of a modification of gravity to- 
day. If the scalar <f> is taken to be relevant only dur- 
ing inflation, however, the Ci are unrestricted and must 
be determined independently from cosmological probes. 
Cosmological constraints on C2, C3 and C4 are quoted by 
Ali et al. [3. 



d x yf^g 



|(V0) 2 + ||n^(v^) 2 - ||(v^) 2 {(n0) 2 - (v„v^)(vv^) - \r(v^) 2 } 



+ |§(V</>) 2 {(D0) 3 - 3(D0)(V Al V^)(V' i V» + 2(V / ,V„0)(VV < V)(VaW) - 6G F V"V^V^V^} 



(10) 
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A. Infrared completion of the Galileon 

Models of this type are unusual, because they must be 
viewed as effective field theories both in the usual sense 
of being valid below some energy scale A, and also in the 
sense that they may require a non-trivial infrared com- 
pletion. For instance, as described in they typically 
arise as intermediate-scale effective theories correspond- 
ing to massive gravities, and are therefore modified at 
an infrared scale set by the Compton wavelength of the 
graviton. Such models were discussed in Refs. @,|8[. The 
consistency of our present analysis only requires that the 
infrared cutoff is much larger than the Hubble scale dur- 
ing inflation. However, in a full cosmological model ac- 
counting for the subsequent expansion of the universe, it 
is likely that this scale must be many orders of magnitude 
larger. We do not address this issue in this paper. 



B. Inflation in the de Sitter decoupling limit 

To get a sense of the background solutions and new fea- 
tures in Galileon models, consider the background field 
evolution in a de Sitter decoupling limit, defined by the 
limit Mp — > oo where H is kept fixed. (We caution that 
this should not be confused with the decoupling limit 
considered in the next section, which allows for a more 
general FKW background) . This limit is applicable if the 
variation AV in the inflationary potential over the du- 
ration of inflation satisfies |AV/V| < 1. In this limit 
we have a Galileon model living on the background of de 
Sitter spacetime, with scale factor a(t) — e Ht . 

After integrating by parts to obtain the action in first- 
order form, and cancelling any boundary terms generated 
by this process 43], the action for a homogeneous field 
configuration <p(t) can be written 



So — 



d 4 x a 3 < 



c 2 - 2 2c 3 H . 3 9 Ci H 2 



14 



6c 5 iJ 3 



(11) 



According to the discussion of Eqs. dU)-©, the current 
J* can be written 

A 3 

J 4 = (f> (c 2 + 6c 3 Z + 18c 4 Z 2 + 30c 5 Z 3 ) = — . (12) 
where Z is the dimensionless combination 

ry _ H(j> 



A 3 



(13) 



We shall see that Z plays the role of coupling constant 
in Galileon theories. 

There are two regimes. First, a weakly coupled solu- 
tion for which Z -C 1, 



3c 2 H 



(14) 



In this regime the outcome is very close to canonical slow- 
roll inflation. Second, there is a strongly coupled regime 
for which Z ^> 1. Specializing to a DGP-like Galileon 
theory [20( for which only C2 and C3 are nonzero, the 
strongly coupled regime gives 



' A 3 A 3 
18c 3 H 2 ' 



(15) 




The smooth solution interpolating between Eqs. (fT4")) 
and US} is 



(16) 



More generally, when the other Galilean interactions are 
present, the field configuration interpolates between the 
weak coupling regime ([T4"f and the strongly coupled solu- 
tion 4>H ~ (AA 2 ) if c 4 # or $H ~ (A 3 A 9 ) 1/4 if C5 ^ . 

What is the relevance of corrections to (fTTj) from non- 
minimal coupling to the geometry, which were briefly 
discussed below Eq. (fTU)) ? To be concrete we consider 
a coupling of the form (f>G, where G is the Gauss-Bonnet 
invariant, which would retain the important feature of 
second-order equations of motion. Since G involves the 
square of contractions of the Riemann tensor, the lead- 
ing contribution must be of the form ~ H 3 (f>/A. This will 
be of comparable importance to the C3-dependent term 
- H<j) 3 /A 3 only if 



JL£>i 

Z 2 A 4 ~ ' 



(17) 



Unless Z < (H/A) 2 <C 1, Eq. <jT7J) and similar conditions 
for the C4 and C5-terms show that the Gauss-Bonnet in- 
variant is negligible compared to the terms proportional 
to Cj, arising from Galileon operators. On the other hand, 
the relative importance of the higher-derivative Galileon 
operators to the lower-derivative operators is controlled 
by positive powers of Z. 

This has a simple interpretation. When Z > 1 the 
nonlinearities of the Galileon sector are important, as in 
Eqs. (|12j) and (|15[) . In this limit the mixing with grav- 
ity can be neglected, and non-minimal couplings such as 
4>G are irrelevant. It is in this regime that an interesting 
Vainshtcin effect can emerge. In addition, as we will ex- 
plain below, the inflationary phenomenology is rather dif- 
ferent to canonical slow-roll models. On the other hand, 
the limit Z <C 1, as in Eq. (fT4")) . describes weakly coupled 
perturbations in the Galileon sector, giving a theory al- 
most equivalent to canonical slow-roll inflation. In this 
limit the strongest interactions come from mixing with 
gravity, and a coupling such as <pG cannot be neglected. 



C. Effective field theory for inflation 

How docs an era of Galileon inflation differ from a 
canonical inflationary phase? We have argued in Sjlllthat 
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two Lagrangians are inequivalent only if they make differ- 
ent predictions for observables. Therefore we must study 
perturbations generated by the action ([2]), which are the 
appropriate measurable quantities. 

Unitary gauge action. Cheung et al. [3l| argued that 
certain properties of the perturbations generated in an 
inflationary model were fixed by the background, and 
were therefore model independent, whereas others varied 
between theories and could be used to probe different 
choices of Cm- This conclusion was obtained by con- 
structing the most general action for small fluctuations 
on a quasi-de Sitter background, subject to the condition 
of unbroken spatial diffcomorphisms and nonlinearly real- 
ized Lorentz invariance. In i]IVDI below we will construct 
the action for small fluctuations in Galileon inflation us- 
ing a more direct approach. However, this action could 
equally have been obtained by specializing the result of 
Ref. [31] to a scenario with Galilean symmetries. There- 
fore, before proceeding with a detailed calculation, it is 
of interest to determine what constraints are placed on 
the model by the construction of Cheung et al. 



The authors of Ref. [3l| worked in a model with a sin- 
gle scalar field, <fi, and constructed their action in a gauge 
where slices of constant time coincided with slices of uni- 
form <p. In this gauge there are no explicit scalar fluc- 
tuations, but only perturbations of the metric. The unit 
vector normal to slices of constant time is n M , and consti- 
tutes a preferred vector field breaking manifest Lorentz 
invariance. The Lagrangian for inflationary perturba- 
tions can be built only out of operators which are in- 
variant under spatial diffeomorphisims associated with 
reparametrizations of the induced three-dimensional spa- 
tial metric h^ u = g^ v + n^n u . Cheung et al. showed that 
it was sufficient to take the Lagrangian to comprise a gen- 
eral scalar combination of the Riemann tensor, R^upat to- 
gether with the time-time component of the metric, g 00 , 
and the extrinsic curvature, K^ v = — h(,, p h,A C '^ „n„, as- 
sociated with slices of constant time [31(. Bartolo et al. 
argued that the most general Lagrangian including terms 
of up to cubic order in small fluctuations can be written 

M 



S= d 4 X y/=i 



1 



-M$R - c(t)g m - A(t) + -M 2 {t)\g m + l) 2 + -M 3 (t) 4 (<7 00 + 1) 



Mi(tf 
2 

M 4 (t) 3 
2 

M 7 (t) 



{g m + l)5K» 
(g 00 + l) 2 5K% - 



M 2 {tf 
2 

M b {tf 



3 



5K^5K, LU 



(g 00 + l)(5K\f - M^)l(g°° + l)6K^6K t 



(18) 



where we have used the Riemann tensor only in the form 
of the Ricci scalar, to match low-energy gravitational ex- 
periments which probe the Einstein action. 

Following Refs. [HI [33j], in writing Eq. (1T%1) we have 
organized the expansion in powers of 6g oa = g 00 + 1 
and 5K pv , which are fluctuations around an unperturbed 
Friedmann-Robertson- Walker geometry. Therefore only 
the operators multiplying c(t) and A(t) are non-zero on 
the background, which fixes these coefficients in terms of 
the expansion history H(t). The Mi(t) and Mi(t) are not 
fixed by the background evolution and encode differences 
between models determined by our choice of Cm ■ It is in 
the effects generated by these operators that we should 
look for distinctive signatures of Galileon inflation. 

To convert Eq. (fT5]l into a form suitable for com- 
putation it is helpful to make a gauge transformation, 
t — > t = t — 7r(x, t). After this transformation, the equal- 
time hypersurfaces of the uniform-field ('unitary') gauge 
are deformed by 7r, which Cheung et al. argued could be 
considered as the Goldstone boson associated with bro- 
ken time translation invariance. Each term in Eq. (|18[) 



generates a series in powers of n, with each copy of tt 
carrying least one gradient in time or space. Each copy 
of 5K pu may additionally contribute a single power of 7r 
carrying two gradients. It follows that Eq. (|18p gener- 
ates cubic interactions involving three copies of 7r with 
between three and six derivatives. In this way a descrip- 
tion of the system is constructed in terms of the low- 
est dimension operators compatible with the underlying 
symmetries — the effective field theory approach. 

Decoupling limit. Eq. (fT5|) is quite generally applicable 
and accounts for the mixing between 7r fluctuations and 
the metric. However, when large nonlinearities are asso- 
ciated with the self-interactions of 7r there exists a decou- 
pling limit in which reliable predictions can be extracted 
while neglecting gravity 0, [2?], Hi], [3^] ■ For this limit to 
be a reasonable approximation, the Mp-suppressed terms 
which are neglected must be smaller than terms which do 
not vanish when Mp — > oo. 

We briefly recapitulate the argument of Ref. [3l| • The 
most relevant kinetic term for the metric fluctuation Sg 00 
will come from the Ricci scalar. Therefore we can pass to 
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canonical normalization by the rescaling Sg 00 —> 5g®° = 
MpSg 00 . The most relevant kinetic term for ir will arise 
from some nonlinear operator in Eq. (|18)1 . If it is minimal 
in derivatives it will be of the form M 4 tt 2 , where M is 
some combination of the scales Mi or Mi. The canoni- 
cally normalized field is tt c — M 2 ir. At quadratic level, a 
mixing term such as M i 7rSg 00 is negligible in comparison 
with the 7r kinetic term for wavenumbers k which satisfy 
k > E mix = M 2 /Mp. The same applies for cubic terms, 
where the leading mixing term M 4 Tr 2 Sg 00 is negligible 
in comparison with M 4 tt 3 under the same condition. A 
similar argument can be given if the most relevant ki- 
netic term for n contains higher derivatives [3l| . or if 
the leading cubic terms enter with a mass scale different 
from M. In the decoupling limit the scale -E m i X must 
be smaller than the Hubble scale during inflation, mak- 
ing our predictions accurate to a relative error of order 
E m i K / H . In what follows we work in this limit, in which 
the metric can be taken to be unperturbed. Therefore 
it is most convenient to work in the uniform curvature 
gauge, where the unperturbed metric is spatially flat and 
can be taken as the background de Sitter geometry. 



the time-dependent coefficients a, f3 and /; satisfy 



D. Fluctuations 



On the basis of the foregoing discussion, we should 
study the effect of fluctuations around a cosmological 
background solution in the decoupling limit by construct- 
ing small fluctuations t M> t + £(x,t) on a hypersurface 
of constant time, working in a gauge where such hyper- 
surfaces are spatially flat. Note that this gives £ an en- 
gineering dimension [mass] -1 . 

The comoving curvature perturbation, £, satisfies £ = 
Hi; and will be conserved on superhorizon scales, where 
kc s /aH — > 0. Working to cubic order in £, the action can 
be written 



SD 



ft ^ ft ^ n ^ 



(19) 



where the symbol 'Z>' is used to denote that S contains 
these contribution among other higher-order ones, and 



a = y (c 2 + 12c 3 Z + 54c 4 Z 2 
/?= y<jc 2 +4 C3 (2Z+^) 



120c 5 Z 3 



(20) 



A 
h 
h 




(21) 

(22) 
(23) 

(24) 



30 3 



A 6 



-( Ci + Ac 5 Z). 



(25) 



The quantity Z was defined in Eq. (1131) . In order that 
Galileon self-couplings dominate the interactions, and 
mixing with gravity can be neglected, we must have have 
Z > 1. Although one can contemplate the limit Z ^> 1, 
there is some risk that this would spoil inflation unless 
renormalized by a Vainshtein effect. In this paper we re- 
strict our attention to the case Z ~ 1 where nonlinearities 
are significant but not problematic. 

The contribution proportional to £ 2 <9 2 £ can be removed 
after a field redefinition. Making the transformation £ — > 
it = £ + f 2 i 2 /2/3, it follows that Eq. (JTHJ can be written 



SD d 4 x a 3 



a < 7r 



gin 



(26) 



where g s = / 3 , g 4 
gi is defined by 



fi, we have defined c 2 = /3/a, and 
13 d ff 2 



2H 2 a 
9i — Ji H + rr-/2 - „ 2 , 

ci Sci a ici at 



13 



(27) 



Had we worked from the uniform-field gauge action of 
Cheung et at, Eq. (|18[) . we would have obtained a, c s 
and the gi in terms of c(t), A(i) and the theory-dependent 
scales Mi(t) and M{t). 

Although £ and 7r will differ on small scales, they be- 
come equal whenever £ = and are therefore equal in any 
epoch when £ is conserved. In particular, they coincide 
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on superhorizon scales. Therefore, to obtain the corre- 
lation functions of £ it suffices to obtain the correlation 
functions of 7T. 

Relation to EFT action. The Galileon Lagrangian, 
constructed above, could have been recovered from 
Eq. (|18p by imposition of the Galileon symmetry, after 
taking advantage of possible field redefinitions and inte- 
gration by parts. Therefore, one might have some reser- 
vations that Eq. (p?S|) is in conflict with the conclusions 
of Cheung et al, who found that the coefficient of the 
operator ^(dn) 2 was fixed by Mi(t). In a generic model 
this coefficient is also responsible for a nontrivial speed 
of sound, c s < 1. Therefore the coefficient of 7t(<97t) 2 is 
fixed once the background evolution and c s have been 
specified. On the other hand, Eqs. ([2l))) - ([2"Tl) show that 
the coefficient of -h(dn) 2 in the Galileon theory is inde- 
pendent. 

This apparent discrepancy disappears if one accounts 
for all terms in Eq. (I18p . which in principle contains 
fourteen free coefficient functions. Of these, the Planck 
mass is measured by terrestrial experiments and the pair 
{c(i), A(i)} must be chosen to match the expansion his- 
tory H(t), leaving eleven free coefficients overall. Con- 
tributions to 7r((?7r) 2 arise from many of these operators, 
which are made relevant owing to the symmetries of the 
Galileon theory. (See £HVE[ where the relative magni- 
tude of each term is clearly expressed by their contri- 
butions to the observable parameter /nl-) It is these 
additional terms which break the expected correlation 
between g% and c 2 . 



E. Primordial density perturbations 

We now proceed to compute the form of the bispec- 
trum for Galileon inflation. As we have explained above, 
we focus on nongaussianities because we expect them to 
parametrize the difference between inequivalent choices 
of the inflationary Lagrangian. Although it is also im- 
portant to study the properties of two-point statistics, 
these are effectively constrained to match observation by 
the closeness of the background solution to a de Sitter 
era with slowly varying H . 

When computed using Eq. (|26p . the correlation prop- 
erties of the primordial density perturbation can be ex- 
pressed in terms of the coefficients of the relevant oper- 
ators, which are a, c s , gi, g% and 54. These represent 
the largest contribution to each observable. We have ar- 
gued that there are two irrelevant operators, X 3 cf> and 
m 2 </> 2 , which although not invariant under the Galilean 
symmetry can self-consistently be made small, but play 
an essential role in ending inflation. These correct the 
predictions of Eq. (|26[) . Their influence can be accom- 
modated by inclusion of the first subleading slow-roll cor- 
rections. As remarked by Kobayashi et al. [1^ . another 
reason to calculate these corrections is that c s can exhibit 
a modest but non-negligible variation over the duration 
of inflation, which manifests itself as a correction which 



is formally of first subleading order. 

Corrections of this kind were calculated for the two- 
point function in canonical inflation by Stewart & Lyth 
[44|. The two-point function describing correlations in- 
duced by the quadratic part of Eq. ([25]) was given at 
leading order by Garriga & Mukhanov 27]. Sublead- 
ing slow-roll corrections were later obtained for the case 
£m = P(X, 4>) by Chen et al. (l8j |. who also gave expres- 
sions for slow-roll corrections to the bispectrum. These 
involved quadratures of the sine and cosine integrals Si x 
and Cix which could not be evaluated in closed form. 
Very recently, Kobayashi et al. |l2| obtained slow-roll 
corrections to the two-point function in models of the 
form Cm = P{X, 4>) + F(X, </>)□</>, which includes the 
term involving C3 in Eq. (110[) as a special case but not 
the terms containing C4 and C5. 

Two-point correlations. We define s = Hc s /c s [24| . 
which measures the time dependence of the speed of 
sound. Similarly we shall require the rate of variation 
per e-fold of each time-dependent coefficient in Eq. (|2l)j) . 
We define 



v = 



Ha 



H 9i 



(28) 
(29) 



for i = 1, 3 and 4, and treat all these combinations as 
the same order of magnitude as e and r\. Certain combi- 
nations of these parameters occur frequently, for which 
it is convenient to define abbreviations, 



A = e + 



v 3s 



^0 = e + v + 2s + i— A 
Mi = e + s ~ 



(30) 

(31) 
(32) 



We note that the combination ^0 + Ml = 2A is purely 
real. 

Translating to conformal time, defined by r = 
J dt/a(t), the two-point function can be written 

(7r(k ls r)7r(k 2 , r')> = (2^) 3 5(k! + k 2 )G fel (r, r') (33) 
in which Gk is defined by 



r - J U t(T')Uk(T) T<T' 

Gk -\ul(T)u k (T>) r'<r 



(34) 



The elementary wavefunction Uk out of which this two- 
point function is built satisfies 

where Hi, is the Hankel function of the second kind. 
The order can be written v = 3/2 + A (45j . 
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The power spectrum, P(k, r), is defined by the equal- 
time correlation function, 



P(k,T)=G k (T,T). 



(36) 



On superhorizon scales, where \kr\ <C 1, Eqs. (j3"3"1) - 
([31)| imply that the power spectrum achieves a time- 
independent value 



P(k,T) P(k) 



H 2 1 + 2E k 

2(kc s \ k ) 3 2a k 



(37) 



where E is the combination 



E : 



-e-s+(2- 7E -ln2)A, 



(38) 



in which 7e is the Euler-Mascheroni constant, and a 
subscript k denotes evaluation at the time kc s r = — 1. 
Eq. (|3"T|) has engineering dimension [mass] -5 , showing 
that the two-point function, Eq. (|33l) , has engineering di- 
mension [mass] -8 which is correct given our definition of 
7T. The time-independence of P(k) outside the horizon 
can be understood as a consequence of the shift symme- 
try of Eq. ([2|), whi ch g uaranteed the conservation of £ 
outside the horizon |4q ]. 

In writing expressions valid to subleading order in slow- 
roll parameters we must be cautious when specifying the 
time at which background quantities such as H, c s and 
the slow-roll parameters are to be evaluated. A small 
shift in the evaluation point translates to a difference 
of 0(1) in the slow-roll suppressed terms. We therefore 



have some freedom to rearrange coefficients by judicious 
choice of the time of evaluation. 

Wavefunction corrections. At the level of the three- 
point function, slow-roll corrections arise from several 
sources. These were identified in Ref. [llj]. There are 
corrections due to the presence of time-dependent, slow- 
roll suppressed factors at each vertex. In addition there 
are corrections involving the combination E which arise 
from taking the \kr\ <C 1 limit of each external line. A 
third class of corrections arise from modifications to the 
wavefunctions carried by internal lines. These were ob- 
tained explicitly by Chen et al. [18j, but in this section 
we give a slightly different treatment which will enable us 
to evaluate the final three-point function in closed form. 
Further details are given in Appendix IA1 

The background wavefunctions can be written 



u(k, t) 



iH k 



(fcc s | fc ) 3 / ; 



■(1 - ifcc s | fe r)e 



i/,r. 



(39) 



The 0(e) correction, 8u{k,r), is obtained by expand- 
ing Eq. (|35[) uniformly to first order in small quantities. 
The variation with respect to the order, is, of each Han- 
kel function can be evaluated using expressions (B.42)- 
(B.46) of Ref. [l8j]. Assembling sine and cosine integral 
terms, these expressions can be rewritten to find 



dH [2 \x) 



i 



dv x 3 / 2 
e i:E (l-ia;)Ei(-2ia; 



2e" 



ife-- ( l 



(40) 



Using Eq. (|4"0l and Eq. (I3"5l , and rotating the contour 
of integration of Ei, we finally obtain 



Su(k, r) 



igfc 1 

(*c s | fc ) 3/2 



\ k e- l ^ T (l + ikc 3lk r) I ^e 2i 



(41) 



jkc s t k T 



Mo|fc + Wi\kkc s \ k T + s k k 2 c 2 s \ k T 2 + X k N k - i\ k kc s \ k N k T - s k k 2 c 2 s ^ k NT 2 



We have defined N k = ln|/cc s | fe r|. The remaining inte- 
gral is to be taken over a contour displaced slightly above 
the negative real axis for large |£|, which renders it fi- 
nite. Similar integrals are generated in the Schwinger (or 
"in— in") formulation of quantum field theory, where the 
same contour prescription is obtained after accounting 
for ie terms which project onto the vacuum at past infin- 
ity [29l . l47j . Differentiating with respect to r and using 
Eqs. (13111) - ([52"|) one can confirm that, despite the appear- 
ance of an apparent logarithmic singularity, Su(k, r)' — > 
in the limit r — > 0. This is the same behaviour as u(k, r) 
itself, and guarantees that the introduction of slow-roll 
corrections does not cause a convergent time integral to 



become divergent. 

Log-divergent integrals of the form appea ring in 
Eq. (|4T1) have previously been obtained in Refs. [3fll48l. 
|49[, which discussed the possibility of singularities for 
certain kinematic configurations of the fcj, including the 
'squeezed' configurations where one ki becomes much 
smaller than the other two. The behaviour of the three- 
point function in this limit is not trivial, but it can be de- 
termined nonperturbatively in the single-field framework 
following an argument due to Maldacena [2{| [5(| , and is 
known to be regular. Therefore any singularities arising 
from this log-divergent integral must cancel. We have 
confirmed that our final expressions contain no singular- 
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ities, but we discuss the significance of these potential 
divergences in Appendix IA"1 

Three-point correlations. We give technical details of 
the calculation of the three-point functions arising from 
each cubic operator in Eq. (pZfH) in Appendix [B] In this 
section we report the final values of /nl, specialized to 
the equilateral limit where all fcj have a common magni- 
tude. We would typically expect the bispectrum to be 
maximized on a configuration close to equilateral, and 
this limit should give a good estimate of the magnitude 
of the bispectrum on this peak configuration. 

We adopt the convention that background quantities 
are to be evaluated at the horizon-crossing time corre- 
sponding to the symmetric point kt = k\ + fc 2 + k 3 , and 
denote evaluation at this time by a subscript This is 
somewhat larger than any individual fcj. In the equilat- 
eral case this moves the point of evaluation to In 3 ~ 1 
e- folds after the common time of horizon exit. 

We define the bispectrum B T (k\,k 2 -,k 3 ) by 



make B a function of the magnitudes fcj alone, indepen- 
dent of the relative orientation among the k^. We define 
/nl to be the reduced bispectrum 



B(fci,fc 2 ,fc 3 ) = |/nl 
5 

x P(h)P(k 2 ) + P(h)P(k 3 ) + P(k 2 )P(k 3 )] , 



(43) 



where all quantities are evaluated at time r. Our con- 
vention that the background quantities in each copy of 
the power spectrum P(k) are to be evaluated at r* im- 
plies that there is a logarithmic correction to Eq. (|3"T)l 
proportional to \nk/k t . Combining Eqs. (|37|) and (j43]) 
we obtain 



Ilikf 



6, = 



B(k 1 ,k 2 ,k 3 ). 



4a*c;j 



(44) 



(7r(ki,T)7r(k 2 ,T)7r(k 3 ,T)) = 

(2^) 3 ( 5(k 1 +k 2 +k 3 )B T (fc 1 ,fc 2 ,fc 3 ) 



(42) 



where the momentum-conservation condition allows us to 



To obtain our final answers we expand this expression 
uniformly to first order in quantities of 0(e). We de- 
fine a numerical constant u>, satisfying cothcj = 5. In 
the equilateral limit fcj = k for all i, each /nl becomes 
independent of k for dimensional reasons and we find 



6 , 3 2 ff i*#* / 2 7E - 3 , 160w - 2 7E - 29 480w - 98 320w - 2 7E - 63 \ 
5 /nl = 27^T I,' + —2~ hl * + 2 + 2 S * + 2 ) 

' (1 - 0.923/11* + 1.141u* - 0.344s* + 0.360e*) (45) 



27 a* 



QMdnf 17 g 3 A , 17 7E -9, 64w-17 7E + 22 192w-34 7E + 40 
5 /nl = - 54cL a, I 1 + -^T~ h3 * + 17 V * + 17 

128w - 17 7E + 18 \ 
+ 17 £ *J 

17 g 3 A 



1 + 0.048/i 3 * + 1.480u* + 3.489s* + 2.008e*) (46) 



54c^ a* 

6 rd^{d.f 13 gt*Hl ( 6-te - 5 173 - 26 7E - 512w 733 - 156 7E - 2304w 

- L = 1 H hA* H v* H 

5 JNL 27ci a* \ 6 26 39 



13 

13 gA*Hl 



147 - 26 7E - 512a; 

e* 



27cl a* 



(1 - 0.256/i 4 * + 2.084v* + 4.509s* + 2.169e*) . (47) 



On the other hand, in the 'squeezed' limit where one fcj 
becomes much less than the other two, each /nl decays 
to zero. Expressions for /nl which describe the complete 
momentum dependence can be extracted from the three- 



point functions given in Appendix [Bj but because they 
are lengthy and ultimately not illuminating we do not 
write them explicitly. 

The leading-order contributions from each of these op- 
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erators were recently computed by Mizuno & Koyama 
fl3j . We have verified that the leading terms of Eqs. (|B1I) , 
(|B9jl and (|B25I) . which generate the leading-order terms 
of Eqs. ([35 ]> -p7 jl . correspond with Eqs. (32)-(36) of 
Mizuno & Koyama. 

Note that, although we are using the conventional no- 
tation '/nl' to denote the reduced bispectrum, Eqs. (|45l) 
(|4"7)) are not directly measurable quantities, because (|42l) - 
(|4"3")l define them in terms of ir. Indeed, according to 
these definitions, Eqs. (|4"5]) -(l4"71) express (ir 3 ) / (ir 2 ) 2 , af- 
ter removal of the momentum-conservation (5-function 
from each correlator, for different choices of (tt 3 ). There- 
fore they possess engineering dimension [mass]. The ob- 
servable quantity is the ratio (C 3 )/(C 2 ) 2 - The appro- 
priate /nl which measure this ratio are obtained from 
Eqs. (|4"5)) ~([4"T |) after division by if*, and are dimension- 
less. 

Dependence on c 2 . Eqs. (j45|) - (j47l) express predictions 
for the three-point correlations generated by the Galileon 
Lagrangian, assuming that the coefficient of each relevant 
operator can be determined by measurement. 

This pattern of three-point correlations gives rise 
to an interesting phenomenology, considerably broader 
than has previously been encountered using noncanoni- 
cal models. In theories such as DBI and fc-inflation, it is 
a familiar result that /nl ~ c^ 2 [18]. Eq. (|47| already 
shows that this is not guaranteed in a model exhibiting 
Galilean invariance, unless the background conspires to 
require 94 /a proportional to c 2 , and as we will argue 
below this is not automatically the case. This effect is 
visible in the calculation of Mizuno & Koyama, although 
these authors did not discuss its significance. 

The possibility that /nl is not proportional to cj 2 in 
the limit of small sound speed has not previously been 
noticed. Why is this? Cheung et al. observed that the 
leading relevant operator contributing to the three-point 
function would be 7r(<97r) 2 , since this is suppressed by 
fewest gradients. In a generic theory where each op- 
erator enters with an approximate common mass scale 
M, the operator 7r(<97r) 2 is dimension six after canoni- 
cal normalization, whereas d 2 ir(dir) 2 is dimension seven. 
Accordingly we would expect 7r(97r) 2 to dominate /nl at 
wavenumbers k < M. The crucial point is that, if ■h(dn) 2 
is the only relevant cubic operator, then it arises from a 
term which also fixes the speed of sound 31]. As a re- 
sult, the nonlinearly realized Lorentz invariance requires 
/nl ~ c~ 2 . 

The important point we wish to emphasize is that 
this is not forced to occur in Eq. (|47|) . because it is a 
dimension-seven operator rather than n(dir) 2 which gives 
this contribution. We will discuss this in more detail be- 
low. One might wonder whether even more powers of c 2 * 
can accumulate in the denominator. However, this is not 
possible. The factor arises from an operator which is 
not suppressed by any temporal gradients, each of which 
contributes a factor c 2 *. No matter how many spatial 
gradients are added, no greater enhancement is possible 
at small c s +. Before drawing any conclusions about the 



scaling of /nl with c s +, however, it is necessary to spec- 
ify the unknown coefficients and a*, and for this one 
must work with a concrete model. 

Mizuno & Koyama discussed two such models, each 
of which could be written in the generic form Cm = 
P(X, <fi) + F(X, </>)□</>, including an approximation to the 
DBI Galileon example introduced in Ref. @ ■ They stud- 
ied this example in two limits, differentiated by &d <C 1 
and bo 1 in their notation, respectively corresponding 
to the dimension-six and dimension-seven operators dom- 
inating the cubic interactions. In each of these extreme 
limits they found /nl ~ cj 2 , because in these limits the 
coefficient of the dominant operator also determines c 2 . 
In the limit where the dimension-six operator is domi- 
nant, the dimension-seven operator gives a contribution 
to /nl of order which is reproduced by our formulas. 
However, on its own this is of limited interest because 
this contribution must be subdominant overall. 

The full, covariant Galileon result extends this in an 
interesting way. First, it shows clearly that it is not nec- 
essary to assume that one operator or the other is dom- 
inant. As we will explain, the Galileon symmetry makes 
it natural for them to be of equal importance. 

Second, the most general Galileon model has three 
higher-derivative interactions, each of which contributes 
to the speed of sound, c s . If any one of the C3, C4 or 
C5 operators dominates a and /?, then c 2 reduces to a 
simple rational fraction which is not especially small. If 
the same operator dominates the gi, then each /nl in 
Eqs. (|4"5l) - (|4"T|) reduces to ~ if, giving an /nl for C of 
order unity. To obtain a very small speed of sound, one 
must suppose that the Cj are arranged in such a way that 
/3 becomes small relative to a. Having done so, in the 
absence of other constraints, there is still sufficient free- 
dom to balance the gi in such a way that the dominant 
contribution to /nl scales parametrically faster than c~ 2 
in the limit of small sound speed. This is one of the key 
results of this paper. 

Relevance of dimension-seven operator. Let us 

explain in more detail how the dimension-six and 
dimension-seven operators are naturally of comparable 
importance. 

First consider a generic theory, without a Galilean 
symmetry. If we tune k ~ M, then a window exists in 
which the contribution from both operators may be com- 
petitive, but in this window we are close to the cutoff of 
the effective theory unless it is forced to higher scales by 
a Vainshtein-like effect [4l|. In practice, if c s is very small 
then the contribution from d 2 ir(dir) 2 could receive an ad- 
ditional enhancement, which might make it competitive 
with 7r(c?7r) 2 . However, c s cannot be made too small with- 
out encountering undesirable stability problems [3l], [5l[ . 
Therefore the dimension-seven operator seems effectively 
irrelevant. 

The situation changes if the Lagrangian exhibits a 
Galilean invariance. In this case, the dimension-seven 
operator gives a significant contribution because the 
dimension-six operators explicitly break the Galilean 
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symmetry. Thus, they arise only because there is a non- 
trivial cosmological background which breaks the symme- 
try, as discussed below Eq. (fT0| . Because this breaking 
is done by the background, it is suppressed by powers of 
H/A. This makes the dimension-six operators formally 
of the same order as the dimension-seven ones. This pos- 
sibility has not been covered in previous discussions of 
the effective held theory of inflationary perturbations. 

To see this in detail, Eqs. (|2"2" ]) -([2l) ]) show that 7r(V7r) 2 
and (9 2 7r(<97r) 2 do not enter with a common mass scale. 
Instead, the coefficient of ^(dir) 2 is suppressed by H, 
which is the same order of magnitude as the extra gra- 
dient carried by d 2 n(dir) 2 . Therefore, the covariant 
Galileon model (specialized to de Sitter space) tunes 
these operators to give precisely competitive contribu- 
tions when c s ~ 1. In terms of observables, the same 
conclusion is easy to obtain from Eqs. (H5j) - (jT7)) . in which 
the dimension-seven contribution from Eq. (|47[) is naively 
suppressed by 0(H) in comparison with the dimension- 
six contributions from Eqs. (H5|) -(P5 ]) . However, the 
Galileon symmetry forces gi* and <?3* to contain an ex- 
tra power of H+ [cf. Eqs. ([212]) and (|2"3|)]. making the 
contribution to /nl from each of Eqs. (l45l) -(|47 |) precisely 
comparable. 

From the point of view of an arbitrary effective field 
theory, such tuning would be highly unnatural. It is re- 
markable that in a Galileon theory this tuning is au- 
tomatic, stable and technically natural: this also fol- 
lows from the dimension-six operators' violation of the 
Galilean symmetry. Therefore we may choose them to 
be small while preserving technical naturalness, in the 
precise sense of 't Hooft |52j |. 



V. CONCLUSIONS 

Successful inflationary models must obey an approxi- 
mate shift symmetry in order to generate sufficient e-folds 
of inflation. This shift symmetry allows us to add any 
scalar constructed from gradients of the inflaton field to 
the inflationary Lagrangian. However, when adding these 
gradient terms we must be wary of a number of pitfalls. 
Adding arbitrary higher derivative operators to the infla- 
tionary Lagrangian can lead to a loss of unitarity, due to 
the appearance of ghost states. The functional form of 
the Lagrangian need not be protected from large renor- 
malizations, and if the Lagrangian requires more input 
parameters than can be measured then the theory loses 
all predictivity. 

In this article we have presented a model, termed 
Galileon inflation, which avoids these pitfalls. Building 
on the success of the Galileon models of dark energy, 
we have constructed an inflationary Lagrangian contain- 
ing noncanonical derivative operators whose form is pro- 
tected by the covariant generalisation of the Galileon 
shift symmetry. It contains a finite number of operators 
and gives rise to second order field equations, implying 
the absence of ghosts. 



We have constructed the action which describes fluc- 
tuations about the Galileon inflationary solution to third 
order in perturbations. In contrast to previous claims 
in the literature we find that none of the coefficients of 
the terms at third order in the inflationary fluctuation 
are fixed by matching to the background evolution and 
the two-point statistics. Therefore Galileon inflation is 
an explicit example of a new class of higher derivative 
inflationary models where the nongaussianity is not con- 
strained to obey /nl ~ 1/c 2 - 
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Appendix A: Integrals of the Ei function 

To obtain expressions for /nl in closed form, we are 
obliged to compute several integrals over the exponential 
integral function Ei, 



Eiz = 




defined for real nonzero x. If x < the integral is man- 
ifestly well-defined; if x > it must be understood as 
a Cauchy principal value. If x = there is a logarith- 
mic singularity which cannot be removed by the principal 
value technique, so the origin must explicitly be excluded 
from the domain of Ei. The exponential integral is re- 
lated to the sine and cosine integrals by analogues of 
Euler's formula, 

177" 

Cia; + iSix = Ei(ia;) + y (A2) 

Gix-iSix = Ei(-ix)-y. (A3) 

In constructing the 0(e) perturbed wavefunctions in 
ijlVDI we must evaluate Ei(2ifcc s r), which can be defined 
by complexifying t in Eq. (|A1|) and interpreting the in- 
tegral over a suitable contour. This contour passes from 
the region where |i| — >■ oo with n/2 < argi < 37r/2 and 
terminates at x. Because r < 0, for x — 2ikc s r the termi- 
nal point lies on the negative imaginary axis. Application 
of Cauchy's theorem allows the contour of integration to 
be rotated, obtaining the result quoted in Eq. (|4ip. 
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Convergent series representation. In calculating 
/nl we require integrals of the form 



dCe 



i(fci + fc 2 -fc3)c<,C 



y e- :? -' (A4) 



where although Jo is a function of k\ , ki and k 3 we have 
adopted the convention that only the momentum which 
occurs asymmetrically is indicated explicitly, in this case 
k 3 , leaving the symmetric dependence on k\ and ki im- 
plicit. It is convenient to define variables ^3 and #3, sat- 
isfying 



03 = 1-^ and 9 3 = ±-^, (A5) 
fef 2§ 3 



in terms of which, after further contour rotations, 
Eq. (|A4[) can be rewritten 



I (k 3 ) 



fl 3 k t c s 



du e 



e 3 u 



dv 



We have taken the limit r —> 0, which incurs an ex- 
ponentially small error if t is sufficiently late that the 
wavenumber kt is a few e- folds outside the horizon [29j . 
In principle there may be an obstruction to carrying 
out the contour rotation if $3 = 0, which occurs in the 
squeezed limit where k\ or ki approaches zero. In this 
limit, the (-integral in Eq. (|A4[) diverges. As discussed 
in Refs. [lo, one should define this limit by analytic 
continuation, first carrying out the integral with #3 > 
and only subsequently squeezing one of the momenta to 
zero. For this reason we need not worry about subtleties 
associated with rotation of either integral to the imagi- 
nary axis. There is also a logarithmic singularity when 
9 3 = 0, causing the interior u-integral to diverge. 

The u-integral has a convergent series representation 
for all 9 3 u > 0, which is guaranteed to be the case in 
virtue of our assumptions and the domain of integration 
for u. This convergent series representation is inherited 

I 



from e 2v , which is an entire function. We find 

v 



c- 2 " = 7E + ln20 3 tt + f; ( 2d3 T - ( A7 ) 

t — t n ■ n,\ 



n ■ n 



Because Eq. (| A7|) is valid for all u > we may substitute 
in Eq. (IA6|) and integrate term-by-term, which yields 



Io(k 3 



^k t c s ^ n J 



(A8) 



Analytic continuation. The sum converges for | ^3 1 ^ 
1/2. Since Iq is defined by an integral it is an analytic 
function of 9 3 on an open neighbourhood of the positive 
real axis, although excluding the origin 9 3 = where 
we have noted that Eq. (|A6j) exhibits a logarithmic sin- 
gularity. Therefore Iq may be determined by analytic 
continuation of Eq. (| A8|) from any open set where the 
sum converges. We conclude that Iq has the compact 
representation 



(A6) 7 (fc 3 ) 



•d 3 k t c s 



In- 



1 + 20, 



29 3 



■& 3 k t c. 



■ln(l-0 3 ). (A9) 



Eq. (|A9I) correctly reproduces the expected pole as $3 — > 
0, associated with a failure of convergence in the (- 
integral of Eq. (|A4|) . It also correctly reproduces the 
logarithmic singularity as 9 3 — > 0, which corresponds to 
the limit i? 3 — > 1. This is also a squeezed limit, occurring 
when k 3 — > 0, and forces ki — k%. Away from these singu- 
lar points, Eq. ()A9|) expresses Io(k 3 ) as a well-defined an- 
alytic function of 9 3 . In practice these singularities cancel 
among themselves in our final expressions for /nl , which 
serves as a consistency check on the calculation. Note 
that the equilateral limit is ^3 = 1/3 or #3 = 1. 

Repeating these steps enables us to find representa- 
tions for integrals analogous to Eq. (|A4[) , with insertions 
of arbitrary polynomials in the (-integral, 



I m {k 3 ) 



We find 



m g i ( fc l + fe 2- 



di 

— 00 £ 



(A10) 



h(k 3 ) 
h(k 3 ) 

I 3 (k 3 ) 

Im{k 3 ) 



1 



{d 3 k t c s ) [ 
i 

1 



■[t?3 + In(l-i>3)] 
^ 3 (2 + z9 3 ) + 21n(l-^3)] 



{^ktCsY 



[7? 3 (6 + 3^3 + 2^)+61n(l-i? 3 )] 



m+1 



(i? 3 fe t c s ) m + 1 [ 



20 3 (m + 1)! F 



1 1 2 + m 
2 2 



29: 



m!{ 7 E + ln20 3 + ^ (m+1) (O)} 



(All) 
(A12) 
(A13) 
(A14) 



where we have retained the convention of denoting de- pendence on the single asymmetric momenta alone. The 



15 



transcendental functions appearing here are the gener- function, ip( m K For the purposes of this paper, we require 
alized hypergeometric function, F, and the polygamma only I m with m ^ 3. 

I 

Appendix B: Three-point functions 



The rules used to obtain correlation functions from a Lagrangian such as (|26p are discussed at several places in the 
literature, to which we refer for calculational details [29j, |4jj, 153J. We calculate the contribution to (it 3 ) induced in 
turn by each operator in Eq. (|26|) . The possible operators are 7T 3 , 7r(<97r) 2 and d 2 TT(dir) 2 . In a model containing more 
than one of these, their contributions add linearly to the total three-point correlation function. 

7T 3 operator. It is convenient to divide the calculation into (a) a part containing the 0(1) contribution and 0(e) 
contributions from external lines and the vertex; and (6) a part containing the 0(e) contributions from internal lines. 
We use the notation of Eqs. ([50)) - (j3"2")) . (|3"5)) . and (|42p and label the bispectrum generated by parts (a) and (b) as 
and B", respectively. For part (a) we find 



o(a) _ 9i* 24 1 



1 + 3£* - A, In + (e* + /h*)7e - ^(e* + 3ft x *) 



(Bl) 



To evaluate (&) it is first convenient to obtain an expression for the integral 
1 



J = 



kt c s * 



d ( e ifc ^«) [ 7o + i 7lCs ^ + l2 cU 2 + i l3 cie + 744e 4 



(B2) 



where ./V*(£) = ln|fc t c s *^|. This is sufficiently general to encompass the r-integration required for each operator, 
not including integrations of the exponential integral function Ei which we treat separately and are discussed in 
Appendix lAl Carrying out the £ integral, we find 



J 



1 



k t c s 



7i + Si 2 72 + 35 2 673 + 11S 3 24 74 + 505 4 

Z. 5 + 17?, + ZA I 



17T 

~2 



<5n- 



Si 



2S_ 2 
.2 



65 3 



245 4 



ju2 1 u3 ' 1,4 

+ f t 



k-fc k^ k-f- k-fc V 2 y V /c^ 

In terms of J and the integral ^2 define by Eq. (|A12j) of Appendix lAl the (6) contribution can be written 

9i* Hi 6 1 



[ - Jj + iA*<4l 2 (fc 3 )] + cc. + (fc 3 fc 2 ->• fci), 



(B3) 



(B4) 



where J3 3 denotes (|B3[) with the assignments 70 = 7 = <5o = Si = 0, and 

(B5) 
(B6) 
(B7) 
(B8) 

The notation '+c.c' denotes addition of the complex conjugate of the preceding term, and k 3 — > k 2 — > ki indicates 
that Eq. (|B4[) is to be symmetrized over the exchanges k 3 f-> k 2 and ^3 f-> ki, yielding a final expression which is 
symmetric between the labels 1, 2 and 3. 

■k(dTr) 2 operator. We divide the calculation into (a) and (b) parts, as above. The 0(1) contribution and 0(e) 
contributions from external lines and the vertex give 



7 


= S = Si = 


72 


= S* - flii. 


73 


= k 3 s± 


S 2 


= A* - 2s* 


S3 


= -k 3 s+. 



(a) gz* Hi 2 fc 3 2 (k 1 -k 2 ) 



k t c s + 



1 + 3£* - A* In 



kik 2 k 3 
kl 



h{ki +k 2 ) +2fcifc 2 



j. 



Tr(97r) 



3(a) 



(B9) 



+ cc + (fc 3 -> fc 2 -> fci), 
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where Jg^f^ is defined by the assignments 73 = 74 — 5 3 = 5 4 = 0, and 

70 = e* 

71 = -e*(fci + fc 2 ) 

72 = -e*fcifc 2 
60 = -e* - /13* 

<5i = (e* + /i 3 *)(fci + fc 2 ) 
£2 = (e* + h 3 +)k 1 k 2 . 

The 0(e) corrections from internal lines contribute 



(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 
(B15) 



B(b) = g 3 * Hi 2 fcKkr-ka) 



iia 3 



where ,/, 



-iA*{ J (fc 3 ) - i(fci + fc 2 )c 5 *Ji(fc 3 ) - k x k 2 c 2 s J 2 (k 3 ) 

+ I {ki) - i(fc 2 - ki)c s Ji(ki) + kik 2 c 2 s J 2 (ki) 
+ I (k 2 ) - i(fcj - k 2 )c si ,h{k 2 ) + hk 2 c 2 s J 2 (k 2 )} + Jj^ 
+ c.c. + (fc 3 -> fc 2 -> fci), 
3A) is defined by the assignments 74 = (5 4 = 0, and 



(B16) 



7o 


= 2^o* - Mi* + s * 


(B17) 


7i 


= s*fc 3 + (fci + fc 2 )(2^ u - ,uo* - «*) 


(B18) 


72 


= fci/c 2 (3/xi* - s*) + s*(fc 2 + k 2 ) + s*fc 3 (fci + fc 2 ) 


(B19) 


73 


= -s+kik 2 k t 


(B20) 


*0 


= 3A* - 2s* 


(B21) 


Si 


= (2s* - 3A*)(fci + k 2 ) - s+k 3 


(B22) 


62 


= fci/c 2 (2s* - 3A*) - s+(k 2 + k\) - s*fc 3 (fci + fc 2 ) 


(B23) 


S3 


= s+faktkt- 


(B24) 



d 2 ir(dn) 2 operator. Applying the same procedure to the final operator, d 2 n(dir) 2 , gives a contribution at 0(1) and 
including O(e) terms from the external lines and vertex: 



9i 



Hi 4 fc§(k!-k 2 ) 

*43 a 3 c 10 ktU . k 3 



(l + 3£* - A* In 



k? ) 



1 + 



3fcifc 2 fc 3 + k t K< 



k? 



2 



1 + 



3K 2 



11 



hk 2 k 3 
kf 



- 7e + 



17T 



2 + 



2X 2 



6fcifc 2 fc 3 ^ I 



(B25) 



+ c.c + (k 3 -> fc 2 -> fci), 



where we have defined X 2 = fcifc 2 + fcifc 3 + fc 2 fc 3 = J2i<j hkj. Finally we also require the (b) contribution, from O(e) 
corrections to each internal line, which yields 



B (b) 



9i* 



Hi 2 fcKkx-k,) 
l :i <>? <-t Uikf 

-iA*|/ (fc 3 ) - i(fci + k 2 - k 3 )c si ,h(k 3 ) + (fc 3 fci + k 3 k 2 - fcifc 2 )c 2 */ 2 (fc 3 ) - ikik 2 k 3 c 3 sir I 3 (k 3 ) 
+ h(h) - \{k 2 + k 3 - fci)c a */i(fci) + {kik 2 + k x k 3 - fc 2 fc 3 )c 2 */ 2 (fci) - \kik 2 k 3 clj 3 {ki) 
+ h(k 2 ) - i(fci + k 3 - k 2 )c si ,h(k 2 ) + (fc 2 fci + k 2 k 3 - kik 3 )c 2 s J 2 {k 2 ) - i/cifc 2 fc 3 c^/ 3 (fc 2 )| 
+ c.c. + (fc 3 -> k 2 -» fci), 



+ j: 



a 2 7r(S7r) 2 



(B26) 
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where J 3 * is defined by the assignments 

70 = 3^o* (B27) 

71 = fctOi* - 2Mo*) (B28) 

72 = **(*? + fca + + ^ 2 (2^i* ~ Mo*) (B29) 

73 = -s* [fef (fa + fa) + + fc 3 ) + fc 3( fc i + fc 2)] - i^fafafa (B30) 

74 = -s+fafafakt (B31) 
^o = 3A* (B32) 

51 = -3A*fc t (B33) 

5 2 = -«*(*? +k 2 2 + kl) - ZKK 2 (B34) 

5 3 = s* + fc 3 ) + k%(h + fa) + fcf (fa + fa)] + ZKfafafa (B35) 
(5 4 = s+kifafakt. (B36) 



Note that, because this interaction is symmetric apart from the arrangement of spatial gradients, Eqs. (|B25[) and (|B26[) 
are in fact symmetric under permutation of the labels 1, 2 and 3 except for the overall factor fc|(ki • k^), which arises 
from the specific gradient combination appearing in d 2 TT(dit) 2 . 
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